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Abstract 



,-^ ' We perform ab initio plane wave supercell density functional calculations 

a' 

O ' on three candidate models of the (3x2) reconstruction of the /3-SiC(001) 



surface. We find that the two-adlayer asymmetric-dimer model (TAADM) is 



^ ' unambiguously favored for all reasonable values of Si chemical potential. We 

P^ ' then use structures derived from the TAADM parent to model the silicon lines 

^) ' that are observed when the (3x2) reconstruction is annealed (the (nx2) series 



of reconstructions), using a density- functional-tight-binding method. We find 
that as we increase n, and so separate the lines, a structural transition occurs 



'^ ' in which the top addimer of the line flattens. We also find that associated 

^ ' with the separation of the lines is a large decrease in the HOMO-LUMO gap, 

and that the HOMO state becomes quasi-one-dimensional. These properties 

"X 

;_i ' are qualititatively and quantitatively different from the electronic properties 
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of the original (3x2) reconstruction. 
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I. INTRODUCTION 

Silicon carbide has long been studied because of its technological potential for electronic 
devicesa'0. While most semiconductor compounds exhibit only one stable phase at room 
temperature, SiC possesses many polytypes. Of these polytypes, cubic or /5-SiC has demon- 
strated its potential for use in high-temperature, high-frequency and high-power electronic 
devices. The lattice parameter of /5-SiC possesses a lattice mismatch of ~20% when com- 
pared with the lattice parameters of silicon or diamond. This has several consequences: first 
that when /3-SiC is grown as a film on a Si substrate, the interface between the two materi- 
als possesses an associated strain field. Second, an unreconstructed Si or C-terminated SiC 
surface is under compressive or tensile stress. This stress is a driving force for significant 
reconstructions on both surfaces. 

Depending on the stoichiometry of the /3-SiC(001) surface, a large variety of such re- 
constructions has been found. For silicon-rich surfaces (containing more silicon than the 
Si-terminated bulk structure), the series of reconstructions (3x2), (5x2), ...(nx2) has been 
observedH'Q, while for less silicon-rich surfaces the c(4x2)&l3 or (2x1) reconstructions are 
observed"!!!!. For C-terminated surfaces, the c(2x2) reconstruction is favouredl!2l£J. 

It has been found that annealing the highly silicon-rich (3x2) surface at 1100 K for 
several minutes leads to the formation of a grating of very long (~1000 A), very straight 
lines at the transition from the (3x2) to c(4x2) surface0tJc3, with the separation between 
the lines of magnitude ~6 A. Additional heat treatment leads to the removal of this grating 
and the formation of widely separated atomic lineg^i!3i!Zl, which we refer to as the (nx2) 
series of reconstructions. These complex {nx2) reconstructions consist of two subunits: the 
structural elements corresponding to the area between the lines (believed to be the c(4x2) 
reconstructioro) , and the structural elements corresponding to the lines themselves. 

In STM images of these linesEJ'Ej the bright regions appear to be ~9 A wide and the 
separation of these bright features along the line appears to be ~ 6 A. The surface lattice 
vector of the /3-SiC (001) surface is 3.08 A. Thus it is beheved that these lines are formed 



from the same structural components as the sihcon-rich (3x2) reconstruction, of which there 
are several models (see figures 1-3): the additional dimer row model (ADRM)I1§I1j, the double 
dimer row model (DDRM)E3'y and the two-adlayer asymmetric-dimer model (TAADM)^. 
The ADRM has an excess Si surface coverage of 1/3 ML (or one addimer per surface unit 
cell) and has an orientation of (2x3) compared to our other candidate models which have an 
orientation we term (3x2), the DDRM possesses a coverage of 2/3 ML (or two addimers per 
surface unit cell), and the TAADM has a coverage of 1 ML (or three addimers per surface 
unit cell). In the interests of completeness we also mention the single dimer row model 
(SDRM) which is a unit cell which is a 90° rotation of the ADRM and which we find has 
a total energy 2.67 eV higher than the ADRM, and so is thermodynamically unfavorable 
compared to the ADRM; and the two-adlayer-asymmetric rotated dimer model (TAARDM) 
which is a unit cell which resembles the TAADM, but with the top addimer rotated so that 
it is parallel to the second level ad-dimers, and which we find has a total energy 0.68 eV 
higher than the TAADM and so is thermodynamically unfavorable when compared with the 
TAADM. We will not consider the SDRM or the TAARDM further. 

In Section II we discuss the various calculational methods used in our simulations, and 
how we can compare the grand canonical potentials of the various models. In Section 
III, using ab initio techniques, we present the structural results that we obtain for the 
(3x2) models, the thermodynamic stability of each model, the electronic properties of and 
how they compare to experiment. In Section IV, we perform calculations on the series 
of (?7,x2) reconstructions modeled by tight-binding, using the thermodynamically favoured 
(3x2) model as the parent, and present the details of a change in the structure of the line 
when n >7. We then present our conclusions in Section V. 

II. METHOD 



A. First-principles total energy calculations 

We performed the calculations using ab initio density functional theory using the projec- 
tor augmented wave method (PAW) to handle the atomic coresa. We use the Perdew-Zunger 
parameterisationo of the Ceperley-Aldera treatment of the uniform electron gas. We solve 
the Kohn-Sham equations using the Car-Parrinello algorithmEB. The main idea of the PAW 
method is to split the all-electron (AE) wavefunction into three parts; 

\^) = \^)+T.i\<P^)-\mP^\^)■ (1) 

i 

|\E') is the AE wave function and |\E') is a pseudo (PS) wavefunction analogous to the wave 
functions of the pseudopotential method, or the envelope functions of the linear method. 
The \(pi) are a set of AE partial waves within the core region, while the {(pi) are a set of 
smooth partial waves which coincide with the corresponding AE partial waves outside the 
core region. \pi) are projector functions localised within the core region which obey the 
relation {pi\4>j) = Sij. 

As there are several models for the (3x2) surface unit cell we modelled all of them, using 
the surface slab technique with periodic boundary conditions. The simulation cells were 
six and ten atomic layers deep (six layers were used for preliminary checks, ten layers for 
results quoted here); the bottom two layers were held in their bulk configuration while the 
other layers were allowed to relax, with the bottom carbon layer passivated with hydrogen 
atoms. In our calculations we sampled four points of the Brillouin zone, corresponding to 
the F-point, the J'-point (in the [TlO] direction in real space), the J-point (in the [110] 
direction in real space) and the /C-point (see figure 4). For our structural calculations we 
used a plane-wave cutoff of 20 Rydbergs, but for our electronic spectra and total energy 
calculations we increased the cutoff to 30 Rydbergs. The vacuum spacing between slabs was 
set to 6.08 A . 

The convergence of the simulation with respect to vacuum spacing and plane wave cutoff 
and A;-point sampling was checked. It was found that increasing the vacuum spacing did not 



affect tlie resulting relaxed structure. The difference between a relaxed cell at a cutoff of 20 
Rydbergs and 30 Rydbergs was negligible, with minor differences in surface dimer length 
(~0.10 A at most) and almost identical ad-dimer length (~ 0.01 A difference). Dimer 
buckling was found to be slightly increased for the higher cutoff compared to the case with 
the lower cutoff (~ 0.06 A increase). We can thus say that our calculations are converged 
with respect to vacuum spacing and plane-wave cutoff. 

B. Tight Binding Calculations 

Tight binding methods (TB)EJ have the advantage of being computationally much less 
demanding than first-principles methods, while still affording relatively high accuracy. The 
disadvantage originates in the fact that the approximations involved are usually of a less 
controllable nature, and hence are more difficult to improve. Nevertheless, the extreme 
simplicity and sometimes surprising accuracy of TB methods make them a frequently used 
tool in computational condensed matter and materials science studies. 

We have made use of the TB model known as Density Functional Tight Binding (DFTB), 
due to Porezag and coworkers^. This method goes beyond conventional TB schemes in 
several ways. First, an atomic-like basis set is explicitly employed. Because of this, the 
model that results is non-orthogonal. In empirical TB methods it is customary to work in 
terms of an underlying basis set (normally assumed to be orthogonal), but only the matrix 
elements of the Hamiltonian are used, and the basis set is never explicitly constructed. 
Secondly, the model is constructed by directly evaluating the Hamiltonian matrix elements 
within the framework of Density Functional Theory, albeit using a two-center approximation 
(i.e. neglecting environment effects), while in empirical TB models the matrix elements are 
adjusted either to empirical or theoretical data. For a detailed description of the model 
construction procedure, the reader should consult Ref.E3. 

In DFTB, the total energy of the system under study is calculated as follows. The 
relative positions of the atoms in the system determine the values of the Hamiltonian and 



overlap matrix elements. From these, the Schrodinger equation for the electronic problem 
can be obtained in matrix form, and solved by the usual techniques (for example, by matrix 
diagonalization) : 

/ . C'j/J {Hia,jp — EnSiaJl3) = 0. (2) 

Here indices i and j label atoms, Greek indices label basis set functions, and n labels eigen- 
states and eigen-values. Once the eigen-value problem (|^) is solved, the band-structure 
energy is calculated: 

occ 

Ebs = 2 2_, -E^n , (3) 

n 

where the sum extends over the occupied eigen-states, and the factor of two accounts for 
the degeneracy of spin. To complete the total energy of the system a repulsive potential 
contribution is added, 

^rep = E^(|R.-I^*l)' (4) 

where Rj is the position of atom i. This potential accounts for the core-core repulsion and 
the double counting of the electron-electron energy implicit in the band-structure term (see 
Ref.il). 

We have used a DFTB parametrisation consisting of four basis functions {S, P^, Py, Pz) 
per atom, corresponding to the four valence orbitals of carbon and silicon. Similar models 
have been used with great success by Porezag and coworkersE^ to study carbon clusters and 
hydrocarbons, as well as bulk crystalline phases of carbon, amorphous carbonE^, Si and SiH 
clusterg^a, crystalline BNE3, and SiC structuresE^. 

C. Thermodynamical considerations 

We can use the total energies calculated from ab initio calculations to discover which 
of the four candidate models possesses the lowest grand canonical potential. We assume 
that the formation of the different structures is determined by thermodynamic factors but 



we cannot compare the total energies of different models with different stoichiometries. A 
common method of comparing structures is to find the chemical potential usi for the silicon 
adatoms and then compare the grand potentialg23'c3 

Q = E- iVsi/Wsi- (5) 

However, because of the experimental conditions that produce the (3x2) surface and because 
SiC is a two component system, it is nontrivial to calculate //si- Nevertheless we can estimate 
a value for the range of chemical potential that the silicon ad-atoms experience. We define 
Ai^^^jjj^ in terms of /isi by the relation 

//Si = i?sr" + ^E^c'o^, (6) 

where £'g*°™ is the energy of an isolated silicon atom (-3.80769 Hartrees in the PAW for- 
malism); AEq^^ is the corresponding value for the cohesive energy per silicon atom in our 
system. 

We must decide what physical/chemical process determines the value of fi^i (and hence 
AE'^qJj) in our system. We can make several choices. We choose to consider the incorporation 
of Si atoms into bulk SiC. AEq^^ would then become the cohesive energy for a Si atom in 
SiC. However this cannot be calculated directly, because it is impossible to insert a single 
Si into SiC without creating a structural defect of some kind. We can, though, estimate 
the maximum range of /isi that the Si ad-atoms experienceca. We can write the chemical 
potential per unit cell of an ideal bulk system of SiC as 

/isrc' = /isi + ^c, (7) 

where /isi and yUc are the (unknown) contributions from the Si and C atoms. In extremely 
Si-rich conditions, //si approaches the value of bulk Si. We can consider deviations from this 
bulk value A/isi=Aisi — fi^s^^^ (with a corresponding definition for A/ic). The allowed range 
is determined by the heat of formation of the SiC compound 

AHf = /i^^'^ - ^^«"= - ^llf' = (/isi - /i^f ') + (/^c - ^^'c'') = A/isi + A/ic (8) 



The experimental value of Ai^j is 0.72 eVEj. In a Si-rich environment A/zsi = 0, and 
therefore A/zc = Ai^j. Similarly, in a C-rich environment AyUc = and A/zsi = Aifj. This 
means that furthest deviation the chemical potential of Si can make from the value for bulk 
Si is 0.72 eV. We show our results as a function of AEq^^i o'^sr the range from its value for 
bulk Si AE'^'^jj^(Si) to the value obtained by assuming that the cohesive energy of SiC can 
be equally divided between the carbon and silicon atoms, we call this value AE^^^ii^) ■ 

III. 3x2 UNIT CELL 

A. Structural 

It was found from structural calculations that all models we considered were locally 
stable, with silicon adatoms forming addimers. All addimers were found to be asymmetric 
(apart from the case of the TAADM, where the second layer addimers were found to be 
fiat). The structure obtained for the ADRM is in agreement with previous ab m^t^dl§lH3E3, 
and tight-binding workH, with a strongly bound and asymmetric addimer. The structure 
obtained for the TAADM is in agreement with previous workE^. There is an alternating 
arrangement of first layer and second layer addimers as we look along [TlO], with the top layer 
addimer strongly bound and asymmetric, and the second level addimers strongly bound and 
flat. However, the structure obtained for the DDRM is in disagreement with both previous 
ab initial, and tight-bindingM results, as we find that both addimers are buckled and quite 
strongly bound. Our result corresponds to one of the DDRM models studied previously£3 
(the LAFM-DDRM structure). The other DDRM models were all found to be of higher total 
energy than this structure, in contradiction with previous workcHl. The Si atoms of the first 
Si layer of the SiC proper were also found to be dimerised in all unit cells, in agreement with 
previous worka'EOcJ. Table I provides further information about the equilibrium structure 
of each model. 



B. Thermodynamics 

We use calculated ab initio total energies, to find Q from equation (^ between the 
various models. The results presented (see figure 5) are for our most converged simulations 
(30 Rydberg cutoff, 4 ^-points sampled). We first of all start off with the assumption that 
the Si adatoms are in equilibrium with the bulk. As can be seen, all crossover points (that is 
the values of AEq^^ for which Q is equal for two separate models) are outside the maximum 
allowable range of AEq^^i ^^ determined by AHf. We can thus unambiguously apply our 
calculations of Q to determine which model is thermodynamically preferred. Over the entire 
range of allowable chemical potentials, the TAADM model is favored. 

Our results are in agreement with the previous work of Lu et a/.B3, where it is found that 
the TAADM is the thermodynamically preferred model throughout the entire allowable 
range of A£'QQ}j(SiC). Regarding the relative stability of the ADRM and DDRM , we agree 
with other worktZl'El, where it is suggested that the ADRM is the favoured model, although 
the TAADM is not considered. We find that the ADRM is only favored for AEq^^i which 
are close to A£'QQj^(SiC). But there is no evidence for a (2x3) reconstruction which is less 
silicon-rich than the c(4x2) reconstruction. As we find that the c(4x2) reconstruction is only 
preferred over a narrow range of /i when the ADRM is included, (contrary to observations), 
we conclude that the ADRM does not occur in practice, as it would prevent the formation 
of the c(4x2) reconstruction. We can therefore rule out the ADRM. 

We have therefore found from our ab initio calculations of the grand potential that the 
TAADM is preferred over a large range of chemical potential and the other two models can 
be discounted, with the DDRM found to be never preferred and the ADRM only valid for 
a range of chemical potential which is non-physical. 



10 



C. Electronic 

We can calculate the differences (the 'dispersion') between the Kohn-Sham eigenvalues 
calculated at different points of the Brillouin zone, and compare these with recent photoe- 
mission experiments performed on the (3x2) surfaceElcj. Our results, and equivalent results 
in the literature are presented in Table 2. The magnitude of dispersion in the majority of 
cases is approximately equal to the magnitude of dispersion measured^ (~0.2eV). The ex- 
ception is for the DDRM, where we find a large dispersion along T — J. We attribute this to 
our modeling of a perfectly ordered periodic (3x2) surface, whereas the (3x2) surface has 
been observed to be almost perfectly disorderedE3c3, in the sense that there is no observable 
correlation between the addimer tilt in one (3x2) unit cell, and the next (in either the x3 or 
x2 direction). The measured dispersion of the surface state bands is ~0.2 eV along F — J', 
with no measured dispersion along Je3. We find that the only model which reflects this 
anisotropy of dispersion is the ADRM, with the DDRM and the TAADM both possessing 
larger dispersions along T — J and smaller dispersions along F — J', contrary to experiment. 
As seen in Table II, there are some disagreements among the different calculations. We 
especially highlight the discrepancy between our results and those of Lu et alx^ for the ener- 
getically preferred TAADM model. However we find that the HOMO state for the TAADM 
exists mostly in the vacuum of the simulation slab, and suggest that as we use plane waves 
to describe our electron wavefunction, as opposed to the Gaussian orbitals used in other 
worked, our dispersion results more accurately describes the dispersion of the surface state. 

IV. LARGER UNIT CELLS (iVx2) 

We can model the (nx2) series of reconstructions by the mixing of two sets of reconstruc- 
tions, the (3x2) and the c(4x2). We use the thermodynamically favoured TAADM model as 
the parent reconstruction of the lines themselves. We use the favoured MRAD modePea as 
the surface between the linese^. Due to the large size of the reconstructions, ah initio meth- 
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ods are inefficient for modelling these surfaces. We therefore used the non-self-consistent 
density-functional tight-binding method to simulate these surfacesHo. We performed simu- 
lations for n=3,5,7,9 and 11 sampling both the F point and a (221) A;-point mesh generated 
by the Monkhorst-Pack schemeE3. It was found that structural calculations are converged 
for this choice of fc-point sampling scheme. All simulation cells are ten atomic layers thick, 
with the bottom carbon layer terminated with hydrogen to avoid artificial charge transfer. 
The results can be neatly summarised in Figure 6. 

As can be seen, a comparison of ab initio and tight-binding simulations for the (3x2) 
surface unit cell shows that although the addimer bond length ai is different in the two cases, 
with the ab initio addimer length found to be 2.3lA and the tight-binding addimer length 
found to be 2.62A , we find a similar buckling Az, (0.58A and 0.53A respectively). We find 
that as we increase the width of the surface unit cell from n=3, to n=5, the asymmetry 
and length of the weak addimer ai remain approximately the same (this is qualitatively in 
agreement with previous theoretical workl^, wherein it is found that the asymmetry and 
length of the top ad-dimer remains the same). However, we find that for n >7, the top 
addimer becomes symmetric and the addimer bond length becomes shorter with a bond 
length of 2.35A . This matches STM topographs of the lines, which show that the lines are 
composed of symmetric unitsta'E3. The flat addimer of the (7x2) unit cell has a shorter 
bond length than the buckled addimer of the (3x2) or (5x2) because there is a stronger 
a bond and a much stronger vr bond. The bond angles between the top addimer and the 
top-adlayer-to-second-adlayer bonds all become ~ 109.5°, compared to ~ 85° and 125° for 
the buckled addimer case. It was also found that if the addimer is buckled and then left 
to relax in a unit cell where all other atoms are already in their relaxed positions, then the 
addimer remains buckled, albeit with a higher energy than the fiat addimer. That is, the 
structural transition may be kinematically limited. The buckled addimer structure is 0.25 
eV higher in total energy than the fiat addimer structure; however we point out that this is 
less than the average thermal energy of all the surface silicon atoms at room temperature. 
The fiat and buckled addimer structures could both be accessible to the surface. 

12 



Analysis of the bond order (off-diagonal density matrix elements p*) between the top 
adatom and adatoms of the second adlayer show that as the addimer becomes fiat this bond 
becomes much stronger, with larger a and n components. Associated with this stronger bond 
between top adlayer and second adlayer is the loss of an electron from the top addimer. This 
electron is transferred to the surface silicon atoms directly below the addimer as shown in 
Figure 6. 

We also find that coincident with the isolation of the separate lines, there is a change 
in the HOMO-LUMO (highest-occupied-molecular-orbital lowest-unoccupied-molecular- or- 
bital) gap. Analysis of the electronic structure of the unit cell was performed using a (551) 
fc-point Monk;horst-Paclc net. Within our tight-binding formalism, we find that the HOMO- 
LUMO gap of the (3x2) unit cell varies from 1.48 eV at the zone centre, to 1.14 eV at the 
zone boundary. For the (7x2) unit cell, we find that the HOMO-LUMO gap is decreased 
by a considerable amount, from 0.87 eV at the zone centre to 0.30 eV at the zone boundary. 
We also find that the anisotropy of dispersion of the HOMO state changes; for the (3x2) 
unit cell the HOMO state disperses by 0.05 eV along T — J' and by 0.43 eV along F — J, 
while for the (7x2) unit cell we find that the dispersion of the HOMO state is 0.43 eV along 
T — J' and 0.05 eV along T — J . As the HOMO state is associated with the top addimer, and 
as this is now more strongly bound to the second level addimers, this means that there there 
are now strong connections along the line. Therefore the HOMO state of the isolated lines 
shows strong quasi-one-dimensional behaviour when compared with the interacting silicon 
lines that constitute the (3x2) and (5x2) reconstructions. 

We observe three changes that happen together: that the HOMO state becomes quasi- 
one-dimensional in character, that the top adatom of the buckled addimer moves down, and 
that there is a transfer of an electron from the top adlayer to the top layer of the silicon 
carbide crystal proper. We offer one possible rationalisation for this structural transition: 
that as the lines become separated, the HOMO state becomes quasi-one-dimensional. This 
makes the HOMO state interact more strongly along the line, which means that there is an 
extra contribution to the bonds to the addimers of the second layer. This forces the top 
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adatom closer to the surface, and flattens the buckle of the top addimer. As the dangling 
bond state rises in energy as the addimer becomes flat, a state of the top adatom becomes 
depopulated, and an electron leaves the top addimer, to reside in the silicon atoms of the top 
layer of the silicon carbide crystal directly below the addimer. However we are aware that the 
charge transfer which seems to be involved in this structural transition could be inaccurately 
treated, as the tight-binding method we use is a non-charge-self-consistent scheme. 

V. CONCLUSIONS 

We have performed detailed electronic structure calculations on several different models 
of the silicon-rich (3x2) reconstruction of cubic silicon carbide. It was found that the 1 ML 
TAADM model was the preferred model for the (3x2) reconstruction over the entire range 
of allowable chemical potential. Our calculated dispersion values for the TAADM contradict 
experiment however, we find that the only model that matches the photoemission datacfl 
is the ADRM. The mapping between observed surface states and one electron Kohn-Sham 
eigenvalues is not well defined however. We conclude that the TAADM is the theoretically 
preferred model. 

We have also used our DFT electronic structure calculations as the basis for a tight- 
binding analysis of the (nx2) series of reconstructions, that is those reconstructions that 
correspond to the silicon lines observed on the c(4x2) surfaceSlij c3. We find that there is 
a structural transition associated with a critical value of n, where if tt, >7, then the top 
addimer of the TAADM becomes flat. We find that the HOMO-LUMO gap associated with 
the increasing separation of the lines is reduced, so that within the tight-binding formalism 
the {n > 7x2) reconstructions are narrow bandgap semiconductors. We also find that the 
HOMO state associated with the flat addimer is now confined to the line and is strongly 
dispersed along the line, i.e. it displays quasi-one-dimensional behaviour. Our results thus 
show that the electronic and structural properties of the (nx2) series of reconstructions, 
corresponding to widely separated silicon lines, are very different from the electronic and 
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structural properties of the (3x2) surface. Detailed experimental data on the physical 
properties of the (nx2) series of reconstructions is needed to verify these theoretical results. 
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TABLES 
TABLE L Equilibrium bond lengths for the addimers and buckling of addimers (magnitude) 

in our total energy calculations. All distances are in A . See Figures 1-3 for more. 
Model Addimer bond length Buckling of addimer Surface dimer length 

ADRM 2.30 (ai) 0.50 2.59 (di), 2.47 (^2), 2.49 (^3) 

DDRM 2.35 (ai), 2.34 (02) 0.51, 0.53 2.40 (di), 2.41 (^2) 

TAADM 2.31 (oi), 2.44 (02), 2.42 (03) 0.58, 0.01, 0.00 2.41 (di), 2.41 (^2) 

TABLE II. 'Dispersion' of HOMO state along [TlO] (F - J') and [110] (E - J directions, for 

various models. Also shown are values found in the literature. All values are in eV. 

Author ADRM DDRM TAADM 

This work 0.18 {TJ'), 0.00 (TJ) 0.01 {TJ'), 0.34 {TJ) 0.13 {TJ'), 0.21 (TJ) 

Pizzagalh et aS < 0.10 {TJ'), < 0.10 (TJ) 1.00 {TJ'), < 0.10 {TJ) 

Lu et aE 0.20 {TJ'), ~ 0.00 {TJ) ~ 0.00 {TJ'), ~ 0.50 {TJ) 0.37 {TJ'), ~ 0.00 {TJ) 
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FIGURES 
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FIG. 1. Schematic illustration of the DDRM (Anti-ferromagnetic configuration). 
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FIG. 2. Schematic illustration of the ADRM. 



i ^ i 




^Nv^^X^^^^N^^^^ 


[110] 


o#^l -•<>-•-] -©-O 


— *■ 


($r^^^ — Q^r^^^ 


[110] 


Id, dj , 




o*o*o*o*o 


^P Silicon atoms 




(^ Ad-silicon atoms 
(top layer) 




i^Y^ Ad-silicon atoms 
(second layer) 



( ) Carbon atoms 

FIG. 3. Schematic illustration of the TAADM. 
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FIG. 4. Depiction of the surface Brillouin zone and the four special A;-points sampled. 
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FIG. 5. Plot comparing the grand canonical potentials for the three candidate models, the 
ideal (1x1) surface and a (3x2) slice of the c(4x2) reconstruction, as a function of cohesive energy 
per Si atom AE^^^. The zero of the grand canonical potential corresponds to the grand canonical 
potential of the ideal (1x1) surface. Arrows pointing downwards indicates the critical value AEq^^ 
where ^ is equal for a pair of models. Also indicated are the cohesive energies per atom of silicon 
and silicon carbide and the range of AHf, to illustrate the maximum range of allowable cohesive 
energies, for the case where the Si atoms are assumed to be in equilibrium with the bulk. 
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(3x2) ab initio 




(3x2) DFTB 




(7x2) DFTB 



FIG. 6. Illustration of the different addimer structures we obtain for different unit cells and 
different , methods of calculation {ah initio and DFTB). Note that the addimer for the (7x2) is flat, 
compared to the (3x2) unit cell (obtained by both methods). The circle on the (7x2) structure 
indicates where the electron from the top addimer localizes when the transition from buckled to 
flat addimer occurs. 
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